Generation of spin-wave dark solitons with phase engineering 
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We generate experimentally spin-wave envelope dark solitons from rectangular high-frequency 
dark input pulses with externally introduced phase shifts in yttrium-iron garnet magnetic films. We 
observe the generation of both odd and even numbers of magnetic dark solitons when the external 
phase shift varies. The experimental results are in a good qualitative agreement with the theory of 
the dark-soliton generation in magnetic films developed earlier [Phys. Rev. Lett. 82, 2583 (1999)]. 
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I. INTRODUCTION 

Dark solitons have been predicted theoretically and ob- 
served experimentally in different types of dispersive or 
diffractivc nonlinear media, including optical and mag- 
netic systems 0. Recently, dark solitons have been 
generated in Bose-Einstein condensates of dilute atomic 
gases by engineering the phase of the macroscopic 
wave function with a technique known as phase imprint- 
ing, earlier developed to create optical and matter-wave 
vortices . Phase imprinting is a new tool of manipulat- 
ing coherent matter waves, and it is described as shining 
an off-resonance laser on a Bose-Einstein condensate in 
order to create phase steps between its different parts. A 
given phase step defines the parameters of dark solitons 
travelling in the same or opposite directions, and the to- 
tal number of the generated matter-wave dark solitons 
which can be either odd or even 

As a matter of fact, these recent results on the genera- 
tion of dark solitons in Bose-Einstein condensates can be 
compared with and linked to much earlier experimental 
studies of microwave magnetic-envelope spin-wave soli- 
tons in yttrium- iron garnet (YIG) magnetic films Q. 
Indeed, the first observation of dark solitons in YIG 
magnetic films revealed unusual features of the dark- 
soliton generation observed as a change of the total num- 
ber of generated spin- wave dark solitons from even to odd 
with the growth of the input power. These observations 
have later been explained theoretically by employing the 
concept of the so-called induced spatial phase shift 0, 
which is closely related to the concept of phase imprint- 
ing, as was also discussed in Ref. Q- Recently, a further 
attempt to employ the phase manipulation technique for 
generating single and multiple dark magnetic solitons has 
been undertaken in Ref. [3 . 

The purpose of this paper is twofold. First, we ex- 
tend the concept of phase imprinting implemented for 
dark solitons and vortices in Bose-Einstein condensates 
to the field of magnetic solitons, and generate experi- 
mentally single and multiple spin-wave envelope mag- 
netic dark solitons from rectangular dark input high- 
frequency pulses with externally introduced phase shifts 



in YIG magnetic films. Second, we provide a direct ver- 
ification of the theory developed earlier in Ref. . 

The paper is organized as follows. In Section ^ we 
introduce our model which was first suggested in the pi- 
oneering papers and is described by the cubic non- 
linear Schrodinger (NLS) equation for the magnetic field 
envelope. Section lllll summarizes the basic theoretical 
results for the generation of dark solitons by an input 
pulse with a jump across the low-intensity region . The 
main Sec. llVl presents our experimental results which are 
shown to be in a good agreement with the basic theoret- 
ical predictions. 



II. MODEL 

We consider the evolution of a spin wavepacket in the 
form '^{x, t) = u(x, t) exp{i(fcoa;— wqOIj where the slowly 
varying complex envelope u{t,x) is described by the di- 
mensionless nonlinear Schrodinger (NLS) equation Q, 
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where Vg ~ duj/dk is the group velocity evaluated from 
the spin-wave nonlinear dispersion uj(k, |up) at the car- 
rier wavenumber = fc(u;o), D = d^uj/dk'^ is the co- 
efficient of linear dispersion, and = dijj/d\u\'^ is the 
coefficient of nonlinearity. In the derivation of Eq. 
dissipation is neglected 

In the case DN > 0, Eq. ^ has a solution in the form 
of a dark soliton that can be written as follows 131: 



u(a;, i) = Mo tanh^ eyi^{i{Kx — Vlt)}, 



(2) 



where ^ = {t — t^ — x/vs)/tq. Solution ((2J) describes a 
localized dip in the continuous wave (CW) background 
with the half- width tq, the center to, and the velocity 
Vs = Vg + DK . In the phase factor of the dark soli- 
tons, K is the soliton wave number and O = VgK + 
{1/2)DK'^ H- A^juoP is the soliton frequency shift, which 
are the nonlinearity-induced corrections to the wavenum- 
ber fco and carrier frequency ujq. 
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Solution is a special case of a dark soliton with 
a modulation factor A = I, i.e. when the dark soliton 
has a minimum amplitude of u = (the so-called "black 
soliton"). Due to the tanh-function, this soliton has an 
overall phase shift of $ = tt. However, such pulses usually 
do not appear in experimental situations. Instead, a more 
general form of a dark soliton Q should be used, 

u{x,t) = uo{l-A^sech^^)^/^exp{i[kx-Qt+a{^)]}, (3) 

where the phase factor has the form: 



(t(^) = sin 



A tanh ^ 



K = K 



\/T~42 

VsT^A 



-^^sech^C, 



(4) 



9. 



Solution (PI describes a single "grey soliton" with an ar- 
bitrary value of the modulation depth A. For A = ±1, 
Eq. |j2Jl transforms into the "black" soliton 0. 

The important condition for such dark solitons to ex- 
ist is the phase shift given by Eq. that has to be 
present in the carrier wave 0. In the experimental sit- 
uation, there is no such total phase shift and, there- 
fore, only even-numbered symmetric pairs of dark soli- 
tons with equal modulation \A\ < 1 can be generated. In 
these pairs, the phase shifts of the individual dark soli- 
tons have opposite signs and they compensate one an- 
other, so the total phase shift adds up to vanish, = 0. 
Under such conditions, excitation of odd numbers of dark 
solitons seems, in general, impossible 



III. THEORETICAL BACKGROUND 

Recently, it was shown analytically and numerically 
that when an input pulse without any initial phase mod- 
ulation enters a nonlinear dispersive medium, the gen- 
erated localized wave acquires an induced spatial phase 
shift accumulated during its generation . Such a phase 
shift is negligible for large group velocities, e.g., for opti- 
cal solitons in fibers, but it becomes important for spin 
waves in magnetic films. Moreover, if the initial phase 
shift (Fig. nj is included in the inverse scattering trans- 
form 01, the results can be employed to explain the 
specific features of the generation of both odd and even 
numbers of dark solitons observed in experiment Q . 

In order to demonstrate the role of the initial phase 
shift, we summarize the results for the dark-soliton gen- 
eration in the framework of the normalized NLS equation. 



,du d^u 
dt dx^ 
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Equation Q follows from Eq. ^ after rescaling, by as- 
suming the reference frame moving with the group veloc- 
ity Ug, and rcnormalizing the variables t and x. The input 
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FIG. 1: Top: Dark input pulse with amplitude wo, length L, 
and area So (shaded). Bottom: Corresponding variation of 
the externally introduced input phase (not to scale) [loll . 



pulse is characterized by the phase shift (p and the dimen- 
sionless value 5*0, which is defined as 5*0 = \uq\L^{N/D) 
and is usually called the input pulse area as it is propor- 
tional to the carrier wave amplitude uo multiplied by 
(spatial) pulse length L = VgT . The soliton wavenumber 
shift can then be presented as K = 27r/L. 

According to Ref. 0|, the number and symmetry of 
the excited sequence of dark solitons are defined by the 
transcendental equation for the soliton eigenvalue v, 



tan(i^ — (p) 



(6) 



Every solution of this equation for the real eigenvalue i^n 
corresponds to a dark soliton with the modulation depth 
i^n = \/ Sq — v"^/ L, where k„ has the dimension of the 
renormalized amplitude, |mo| \/ I D). The right-hand 
side of Eq. © is always point-symmetric around the ori- 
gin whereas for the tan-function this is only the case when 
= (E>/2 is zero this latter case, Eq. lO yields 

pairs of eigenvalues with the same absolute value but with 
opposite signs. This explains why only even-numbered 
symmetric pairs of dark solitons are predicted by the 
earlier theory |^. Note that a dark soliton with zero 
amplitude minimum is defined by k„ = \uf)\[N / D)^/"^ ^ 
i.e. at z^,i = 0. This is only the case for (f> = -k /2 and odd 
multiples thereof. 

A finite initial value of the phase shift 7^ allows 
to predict the generation of asymmetric or symmetric, 
odd- or even-numbered dark soliton pulses. The typical 
structure of an iV=3-soliton state and the corresponding 
variation of its phase are schematically shown in Fig. |21 
Such a soliton sequence is characterized by the overall 
number of dark solitons and the modulation of each in- 
dividual dark soliton. According to Eq. these values 
arc completely determined by the input pulse area Sa and 
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FIG. 2: Top: Schematic for the generation of A'^=3 dark soli- 
tons. Bottom: Corresponding simplified phase structure. 
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FIG. 3: Number of dark solitons in the plane (Sc^). The 
circle indicates the A^=3-soliton shown in Fig. |5| The sep- 
aration lines between different numbers of solitons are at 
(f) — Mn ±5*0, M = 0,±1,..., the lines of exact symme- 
try are horizontal lines at </) = Nti/2. For even A^, there will 
be pairs of symmetric solitons; for uneven A^, a single soliton 
of 100% modulation depth surrounded by pairs of symmetric 
solitons. 



the total phase shift This suggests the use of the 
parameter plane [Sq,4>) as presented in Fig. |31 Chang- 
ing the pulse length and/or input amplitude is equivalent 
to moving on the horizontal axis. A change in the total 
phase shift is denoted by a shift on the vertical axis. Ev- 
ery point in that parameter plane can be assigned a value 
of 5*0 and 0, therefore it is directly connected to a dark 
soliton of distinct eigenvalues Vn- 

In Rcf . , it was shown that in physical systems with 
low values of the group velocity, an additional phase shift 
is acquired during the process of pulse excitation. How- 
ever, in the general case, the total phase shift in the input 
pulse carrier wave may originate either from the pulse 
forming setup, i.e. it is provided externally, or from an 
intrinsic mechanism in the physical system. In this paper, 
we present a systematic study of the formation of mag- 



netic dark solitons with a variable external phase shift 
(/)ext- Using these data, we verify, in particular, whether 
an acquired phase shift (/)acq is present at all and whether 
it fits the theoretically expected values. 

In essence, the results of Ref. predict that an ad- 
ditional phase shift 0acq acquired during the process of 
dark-soliton generation can be written in the form. 



bT 



(7) 



' D = i'^gT'^ is the (positive) dispersion 
time. The phcnomcnological factor b (introduced in the 
theory and evaluated there to be close to 2) takes into 
account the influence of the finite length of the back- 
ground CW signal. 

Usually, the sign of 0acq would not be important due 
to the symmetry of Eq. ((HJ. However, as we are using a 
phase shifter in the experiment to produce the external 
phase 4>cxt, the relative sign of these two phase shifts is 
important. In the experiments described below, the sign 
of the acquired phase shift 0acq was opposite to that of 
the external phase shift </>ext- Thus, the total phase shift 
of the input pulse is expressed as 



Ototal 



(8) 



where 0acq is determined by Eq. |7J). Note that there 
is always a fixed 'offset' phase shift, even in the limit 
Sq = 0. This means that all dark solitons in magnetic 
systems acquire a phase shift during the process of their 
generation. However, if the group velocity and, therefore, 
Td is large, as we have for the case of dark solitons in 
optics the phase 0acq is small and will remain unim- 
portant. In this case, the dark solitons are always placed 
on a horizontal line at = on the (5*0, 0)-plane which 
is equivalent to the production of even-numbered sym- 
metric pairs of dark solitons. 

If the input pulse parameters arc in a range that 
Td ~ T, the acquired phase shift may manifest itself by 
producing asymmetric and/or odd- numbered dark soli- 
tons, and the points corresponding to these solitons on 
the parameter plane of Fig. 13 will belong to the parabolic 
curve defined by Eq. lO and having a minimum at the 
point Sq=0. This case was realized in thepioneering ex- 
periments performed by Patton's group . 

On the other hand, if the constant and positive exter- 
nal phase shift t^ext > is introduced externally, while 
duration T of the input pulse is kept constant and the 
input wave amplitude (and, therefore. So) is increased, 
the dependence <?!'totai('S'o) ^StS the form of a quadratic 
parabola with a maximum at So=0 [see Eq. (jH])]. As will 
be seen below, it is this case that is realized in our cur- 
rent experiments. Therefore, the points corresponding to 
different output solitonic pulses obtained in our experi- 
ments with a positive external phase shift and increasing 
input wave amplitude belong to parabolas having a max- 
imum at So=Q (see Fig. |S1 below). 
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IV. EXPERIMENTAL RESULTS 

To study the effects produced by an external phase 
shift on the generation of magnetic dark sohtons, we need 
to design an experimental setup with the following char- 
acteristics: (i) the setup should include a nonlinear dis- 
persive waveguide capable of supporting spin-wave dark 
solitons, i.e. the necessary condition for the formation of 
dark solitons, ND > 0, should be fulfilled; (ii) the param- 
eters of the pulses propagating in the waveguide should 
fit the range where the induced phase shift defined by 
the Eq. {Tj) has a large enough value and could be de- 
tected experimentally; (iii) there should be a technical 
possibility to introduce an external phase shift between 
the leading and trailing fronts of the input pulse, simi- 
lar to the phase shift introduced by the method of phase 
imprinting in Bosc-Einstcin condensates or nonlinear op- 
tics Q. 

For our experiments, we have chosen magnetostatic 
spin waves propagating in a thin quasi-one-dimcnsional 
waveguide made of a single-crystal YIG film. Nonlin- 
ear properties of these waves are well understood, and 
the formation of both dark and bright solitons has been 
observed |^. Moreover, the measurement techniques 
incorporating the so-called delay line setup are well es- 
tablished. In our case, it is especially important that one 
can easily introduce an external phase shift in the input 
pulse simply using a microwave phase shifter. The most 
important feature of the magnetostatic waves in the con- 
text of this study is their comparatively low group velocity 
that should result in a substantial value of the induced 
phase shift defined by Eq. ©(see Ref. for details). 

The experiments are performed on a 1.5 mm wide and 
7/^m thick YIG film waveguide which, due to its rel- 
atively small width, could be considered as quasi one- 
dimensional. The film waveguide is fixed in a standard 
delay line setup consisting of two strip-line antennas of 
45/im width and 9.2 mm separation. A microwave fre- 
quency of 5.065 GHz is used. A magnetic bias field of 
1107.5 Oe strength is applied tangentially to the film sur- 
face and perpendicularly to the propagation direction of 
spin waves and, therefore, the conditions for the excita- 
tions of magnetostatic surface waves are fulfilled. The 
value of the magnetic field is chosen so as to place the 
working point frequency in the middle of the spin-wave 
spectrum, letting all Fourier components of the dark- 
pulse signal pass through. 

A variable phase shift in the input pulse is created in 
the following way: a continuous microwave signal pro- 
duced by a sweep-generator is divided into two parts, 
which are then individually controlled by microwave 
modulators connected to the pulse generators. A stan- 
dard microwave phase shifter is introduced in one of the 
modulation channels. Then, the two resulting pulse- 
modulated signals are combined together in a microwave 
mixer device, and are supplied to the input antenna of 
the experimental delay line (see Fig.^J. 

In the experimental scheme described above, the trail- 
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FIG. 4: Schematic diagram of the experimental setup for gen- 
erating dark pulses. The two pulse generators trigger the two 
microwave modulators with a fixed delay. Variation of the 
length of the first pulse provides the dark-pulse duration. The 
phase shifter introduces the external phase shift <^ext. 
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FIG. 5: Timing scheme of the dark-pulse generation setup; ti 
is varied, thus changing the value of the parameter T. 



ing edge of the bright rectangular pulse created in the 
first modulation channel provides the leading edge of the 
input dark pulse, while the leading edge of the bright 
rectangular pulse created in the second modulation chan- 
nel provides the trailing edge of the input dark pulse. 
The phase shift (/)ext is introduced by the variable phase 
shifter inserted in the second modulation channel. Thus, 
the duration T of the formed dark input pulse (see Fig.O 
is determined by the difference of the time delay At be- 
tween the leading edges of both bright rectangular pulses 
and the duration ti of the bright pulse in the first modu- 
lation channel T = At — ti. The duration T of the input 
dark pulses in our experiment could be changed in 1 ns 
steps, and the duration of the CW background before and 
after the input dark pulse is controlled by changing the 
durations ti and t2 in the first and second modulation 
channels, respectively. In our experiments, the values of 
ti and t2 are chosen to be constant and equal, ti = t2 
— 1/is. The input pulse sequences of the form shown in 
Fig. 13 arc produced with a repetition rate of about 100 
kHz. 

The magnitude of the external phase shift introduced 
into the dark input pulse is controlled by the variable 
phase shifter in the second modulation channel. The zero 
point and the scale calibration of the phase shifter are 
established by adjusting it to the points of minimum and 
maximum interference of two overlapping pulses. The 
sign of the phase shift is established by calibrating it at 
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different frequencies. 

The signal is recorded with a digital oscilloscope and 
a transient recorder (500 MHz) through a diode detec- 
tor. The signal received by the output antenna is am- 
plified using a low-noise small-signal amplifier to drive 
the detector diode in its optimum sensitivity range and 
to achieve the highest signal-to-noise ratio possible. The 
output voltage in this range is proportional to the power 
at the input of the detector diode. 

In our experiment the variable parameters are: (i) in- 
put power Pin, (ii) duration of the dark input pulse T, 
and (iii) external phase shift 0oxt- The input power is 
varied in the interval Pin = 1.6 .. . 100 mW, while the 
dark pulse duration is varied in the interval T = 10 ... 30 
ns in steps of 5 ns. We use the following values of the ex- 
ternal phase shift: O.ITtt, 0.427r, O.Gtt, 0.727r, and O.OStt. 

The experimental results are summarized using the pa- 
rameter plane of Fig. 13 The profiles of the output pulses 
corresponding to the fixed values of the input power Pin, 
dark input pulse duration T, and external phase shift 0ext 
are recorded, and the corresponding points are placed on 
the parameter plane Fig.|21 Counting the number of the 
minima in the output pulse profile corresponding to a 
particular set of input parameters, we can determine the 
number of solitons and compare it with the theory. 

However, to compare the experimental results with the 
theory, the dispersion and nonlinearity parameters of the 
magnetostatic waves propagating in the YIG film wave- 
guide should be either measured or calculated. First, 
the group velocity at the working point is determined by 
measuring the time of the pulse propagation between the 
input and output antennae of the setup. This velocity 
is found to be equal to 5.5 • lO^cm/s (±2%) which is 
close to the theoretical estimate of the group veloci ty o f 
5.3 • lO^cm/s (±2%) obtained from Eq. (55) in Ref. 

The value of dispersion D was established by compar- 
ing Vg at several frequencies ±50MHz above and below 
the carrier frequency loq- Using a polynomial fit function 
of Vg{to), the value of D at the working point was deter- 
mined to be -7.2 • lO^cm^/s (±20%). The large error 
is due to the small differences of Vg at small frequency 
intervals. The theoretical estimation of the dispersion 
coefficient D done using Eq. (55) from Ref gives the 
value of —7.1 • 10'^cto^/s±20%. Thus, in the calculations 
below we use the values: Vg = 5.5 • lO^cm/s (±2%) and 
D = -7.2 • lO^cmVs (±20%). 

The value of the nonlinearity coefficient N could not 
be measured directly. Thus, it was calculated using Eq. 
(52) from Ref. O to give = 9.5 • lO^s'^. 

The results of our experiments for the pulse duration 
of T = 20 ns arc presented in Fig. |B1 Two adjustable 
parameters are used to place the experimental points on 
this graph. The first parameter is the coefficient relat- 
ing the input power Pin to the normalized amplitude of 
the spin wave mq, i-e. jwoP = BPi^. It defines the re- 
lation between the input power Pin and the parameter 
Sq used in the theoretical formula {T)) and, therefore, it 
determines the scale of the graph along the horizontal 




FIG. 6: Experimental results for the parameters of dark soli- 
tons with T = 20ns on the plane of 0totai vs. So. The marked 
points correspond to the soliton profiles shown in Fig. |7| 

axis. Similar to Ref. 0, we chose B assuming that at 
the input power of Pin = 175mW, when thermal effects 
begin to manifest themselves, we have a typical value 
of the spin wave amplitude in the film equal to |uo| = 
0.07 (see Ref. ^3)- the conditions of our experi- 

ment this yields B = 27 ■ 10~'^Vl^~^, which is about three 
times larger than the similar value in Ref |1| . The differ- 
ence may be attributed to a better matching between the 
impedances of the supply line and the input micro-strip 
antenna. 

The second adjustable parameter is the phenomenolog- 
ical coefficient b in Eq. {Tj) which describes the effect of a 
finite CW background of the input dark pulse and deter- 
mines the scale of thegraph in Fig. |51 along the vertical 
axis. Similar to Ref. Q, we take b = 2. 

The parabolas presented in Fig. |S1 arc computed using 
Eq. lO, and they correspond to different values of the ex- 
ternal phase shift 0.177r, 0.427r, O.GOtt, 0.727r, and 0.957r, 
so that the lowest parabola corresponds to the smallest 
external phase shift of 0.177r. The values of the initial 
(linear) phase shifts corresponding to the initial point 
Sq — are computed as 

(0total)o = 4>c^t - 2T/Td- (9) 

The points on the parabolas correspond to the points 
where the experimental oscillograms of the output signal 
are recorded. Samples of such experimental oscillograms 
corresponding to the points denoted by letters " a" , " b" , 
"c" , "d" and "e" in Fig.Elare presented in Fig.[7| The left 
row of oscillograms demonstrates that, when the external 
phase shift grows from 0.177r to 0.957r but the input power 
is constant (i.e.. So = 0.757r), the output signal profile 
with (a) two dark solitons transforms into (b) three dark 
solitons, and then again into (c) two dark solitons. 

The right row of oscillograms in Fig. |7| demonstrates 
that with an increase of the input power at a fixed value of 
the external phase shift O.Gtt, the number of dark solitons 
is increasing as the output oscillograms contain cither 



6 



P.^^ 0.75 7t 
it) = 0.90 71 




P^^ 0.75 71 


(b) 


^ =0.55n: 

, 





F|^^ 0.75 n: 


(a) 




R^^ 0.25 7t 


(d) 


^ =0.10 n 






i^^^= 0.55 It 


























FIG. 7: Example of a series of measured dark solitons with 
equal pulse duration of 20 ns. Left row: equal input power 
Pin — — 9dBm (equivalent to So ~ 0.757r), but different exter- 
nal phase shifts (/lext = 0.177r(a), 0.607r(b) and 0.957r(c). Right 
row: equal external phase shift (f>ext = O.GOtt, but different in- 
put power Pin = -ISdBm (d), -9dBm (b) and OdBm (e). The 
soliton profiles correspond to the vertically and horizontally 
aligned points in Fig. |S| 



one, three, or five dark solitons. These results appear to 
be in a good qualitative agreement with the theory . 

We should mention that the number of dark solitons 
observed in the experimental profiles does not always 
match the number predicted by the theory This is 
especially true for very small and very large values of 
the input power. For small values of Pjn, correspond- 



ing to the values 5'o < 0.57r, the dark soliton is not 
properly formed as the time of the signal propagation 
between the antennae is smaller than the so-called " non- 
linear time" Tat = l/(A^|Mp) during which the spin-wave 
nonlincarity could significantly affect the pulse profile. 
On the other hand, for large values of Pin corresponding 
to 5*0 > I.Stt, nonlinear dissipation and other nonlinear 
effects, not taken into account in our simple theoretical 
model, make the dependence of |uoP on the input power 
nonlinear, which prevents us from making a quantita- 
tive comparison between theory and experiment. At the 
same time, it is clear from the experimental data pre- 
sented above that in the range of the intermediate input 
powers, O.Stt < Sq < I.Stt, the theory gives a good qual- 
itative and even reasonable quantitative explanation of 
the experimental data. 

In conclusion, we have extended the concept of phase 
engineering, demonstrated earlier for matter-wave dark 
solitons in Bosc-Einstein condensates, to spin- wave mag- 
netic solitons. We have demonstrated experimentally the 
crucial role played by an externally introduced phase shift 
in the input pulse for the process and outcome of the 
dark-soliton generation. Our experimental results are in 
a good qualitative agreement with the analytical and nu- 
merical predictions of dark-soliton generation, and they 
provide more direct verification of the theory. 
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